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I.  Introduction 

Situations  exist  where  it  is  necessary  to  measure  the  fluence  and 

the  energy  spectra  of  high  energy  neutrons,  En  i  5  MeV.  At  the  present 

time,  the  various  methods  developed  to  do  this  all  have  several  limitations 

associated  with  their  use.  The  highest  energy  threshold,  found  in  the 

commonly  used  threshold  detectors,  is  approximately  3  MeV  for  ^S.  The 

current  state-of-the-art  as  well  as  the  limitations  of  activation  foil 

l  2 

spectrometry  have  recently  been  described  by  McElroy  and  Barra! 1. 

It  is  pointed  out  that:  1)  activation  rate  measurements,  2)  fission 
yield  errors,  3)  evaluated  cross  section  errors,  and  A)  errors  associated 
with  spectral  unfolding  procedure  may  combine  to  produce  disagreements 
of  up  to  30%  in  fission  spectra  measurements  as  compared  with  other 
methods  of  analysis.  The  short  half-life  of  some  of  the  reactions  is 
yet  another  limiting  factor.  Organic  scintillation  spectrometers  are 
able  to  detect  neutrons  with  superior  efficiency  and  their  output  can  be 
unfolded  into  neutron  spectra.  However,  the  spectrometers  are  not  con¬ 
venient  for  use  in  many  situations.  Their  size,  weight,  and  power 
requirements  are  a  limitation.  Also,  there  are  still  problems  in  spectra 
unfolding  techniques,  particularly  at  high  energies. 

One  of  the  more  significant  new  developments  in  the  field  of  neutron 
detection  has  been  made  possible  by  the  discovery  of  a  new  class  of  charged 
particle  detectors  called  nuclear  track  detectors.  The  detecting  principle 
is  based  on  the  delineation  of  individual  paths  of  heavily  ionizing  charged 
particles  such  as  fission  fragments,  recoil  nuclei,  alpha  particles,  and 
low  velocity  protons.  Delineation  is  accomplished  through  a  selective 


chemical  attack  of  the  narrow,  radiation  damaged  region  along  the  particle' 
trajectory.  The  resulting  narrow  channels,  referred  to  as  tracks,  can  be 
enlarged  to  a  convenient  size  for  direct  viewing  with  an  optical  microscope 
Automated  systems  for  track  counting  are  now  commercially  available. 

A  program  is  underway  to  develop  a  high  energy  neutron  detector 

O 

(HEND)  system.  The  development  program  consists  of  (1)  designing 
and  fabricating  a  series  of  neutron  recoil  particle  detectors  from  a 
general  knowledge  of  recoil  interactions,  energetic  ion  transmissions, 
and  track  registration  characteristics  of  nuclear  track  detectors,  (2) 
calculating  the  response  of  the  individual  detectors  to  neutrons  covering 
their  useable  energy  ranges,  (3)  verification  of  the  calculated  detector 
response  by  experimental  measurements  performed  at  several  neutron 
energies,  and  (4)  development  of  a  computer  program  which  will  unfold 
the  HEND  measurements. 

As  mentioned  in  earlier  reports,  the  detectors  are  composed  of  a 
sandwich  arrangement  which  includes  the  following: 

(1)  a  radiator,  consisting  of  a  flat  sheet  of  a  material  containing 
an  element  whose  nuclei  are  elastically  scattered  by  neutrons  and  ejected 
from  the  sheet, 

(2)  a  degrader,  a  thin  film,  usually  of  gold,  which  reduces  the 
energies  of  recoil  nuclei  transmitted  through  it,  and 

(3)  a  nuclear  track  detector,  a  sheet  of  plastic  capable  of  recording 
the  tracks  of  the  energetic  recoil  nuclei  impinging  upon  it.  The 
detectors  have  the  practical  advantage  of  being  small  (-v  5  cm2  total 
area),  light  (several  grams),  stable,  and  independent  of  any  reaction 


half-lives.  The  useful  neutron  energy  range  is  from  -v  5  to  20  MeV.  At 
energies  of  5  MeV  and  below  track  recording  efficiencies  are  low.  Above 
20  Mew  inelastic  interactions  begin. to  compete  significantly  with  the 
elastic  scattering  process. 

II.  Detectors 

The  detectors  are  designed  to  perform  as  a  series  of  threshold  detectors, 
■"he  particular  nertron  threshold  energy  depends  on  the  mass  of  the  recoil 
?.  c’’.”  and  '■he  thickness  of  the  degrader.  For  a  particular  element  and 
; pe,  and  with  .so  degrader,  neutrons  of  some  minimum  energy  will  produce 
tracks  in  the  plastic  with  very  low  efficiency.  As  neutron  energy  increases, 
recoi'  energy  and  range  also  increase.  Therefore,  efficiency,  not  considering 
the  elastic  scattering  cross  section,  continuously  increases  with  energy. 

When  a  degrader  is  introduced  between  the  radiator  and  the  track 
recorder,  recoils  of  low  energy  are  stopped  in  tl.e  degrader  while  those  of 
higher  energy  are  slowed  down.  Thus  a  degrader  of  a  given  thickness  will 
tend  to  shift  the  response  of  the  detector  to  higher  neutron  energies  by 
a  given  amount.  In  this  manner  a  series  of  degrader  thicknesses  produces 
a  series  of  neutron  threshold  energies. 

If  the  radiator  is  thicker  than  the  ranges  of  the  recoil  particles, 
then  the  efficiency  of  the  detector  will  increase  steeply  with  neutron 
energy  until  the  elastic  scattering  cross  section  begins  to  drop  significantly. 
However,  a  thin  radiator  produces  a  detector  with  a  somewhat  different 
energy  response.  When  recoil  particle  ranges  equal  radiator  thickness 
the  response  begins  to  round  off.  Efficiency  falls  considerably  below 
the  thick  radiator  case  for  higher  neutron  energies.  Detectors  of  this 
type  are  useful  in  that  they  sample  lower  neutron  energies  better,  in 
comparison  to  higher  energies,  than  do  the  thick  radiator  types. 
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The  detection  and  accurate  counting  of  recoil  particle  track  densities 
places  stringent  requirements  upon  the  track  recorder  employed.  For  any 
given  recoil  nucleus,  there  is  both  an  upper  and  lower  level  of  track 
registration  sensitivity  which  the  recorder  must  satisfy.  The  lower  level 
is  imposed  by  the  background  problem.  4He  particles  are  produced  copiously 
within  the  radiators  and  plastic  track  recorders  which  have  been  used. 

Some  degrader  materials  tested— notably,  aluminum— also  have  significant 
(n,  4He)  cross  sections.  In  track  recorders  sensitive  to  4He  particles, 
the  track  densities  from  this  source  can  be  greater  than  track  densities 
due  to  recoil  particles.  Since  the  ranges  of  4He  particles  are  also  much 
greater,  the  only  way  these  spurious  tracks  may  be  excluded  is  by  using 
a  track  recorJer  which  is  insensitive  to  4He  particles. 

At  the  same  time,  the  track  recorders  must  be  sensitive  enough  to 
register  recoil  particles  of  the  full  range  of  energies  produced.  Further, 
the  latent  tracks  must  etch  out  such  that  Vy/Vg  »  1,  where  Vy  =  etch  rate 
alonv,  the  particle  path  in  the  track  recorder,  and  V-.  =  bulk  etch  rate  of 
the  track  recording  material.  If  the  above  condition  is  not  met,  two 
significant  things  occur.  First,  track  registration  becomes  a  strong 
function  of  the  angle  of  impingement  of  the  particle  on  the  track  recorder. 
Detector  sensitivity  is  impaired,  angular  dependence  of  the  detector  is 
increased,  and  mathematical  handling  of  the  problem  becomes  difficult. 
Second,  the  tracks  which  are  recorded  may  be  in  poor  microscopic  relief, 
thus  making  track  counting  less  reliable.  This  is  especially  troublesome 
for  those  particles  with  very  short  ranges,  which  includes  a  large  portion 
of  the  recoil  nuclei  produced. 

The  quality  of  the  etched  surface  of  the  track  recorder  must  .Iso  be 


considered.  In  many  materials  etch  pits  form  daring  etching  that  are  not 
caused  by  ionizing  particles.  These  pits  may  be  indistinguishable  from 
short  recoil  tracks.  Where  accurate  track  counting  is  desired,  such 
recorders  are  not  adequate. 

Irradiations  with  14  MeV  neutrons  were  made  to  test  various  detector 

9  1 9  16 

components.  Radiators  of  three  different  recoil  particles  Eh,  C  and  0 
were  used.  Degraders  were  made  of  thin  aluminum  foil  and  gold  films  which 
had  been  vacuum  evaporated  on  to  the  radiators.  The  track  recorders  were 
Lexan  polycarbonate  plastic.  Lexan  is  a  plastic  in  common  use  which  is 

4 

known  to  be  insensitive  to  He  particles  while  recording  the  tracks  of 

12  16 

stopping  ions  such  as  C  and  0. 

The  detectors  employing  Lexan  did  not  yield  recoil  track  density 
measurements  of  the  quality  required  for  this  application.  The  detection 

g 

efficiency  of  the  Lexan  recorder  to  neutron  induced  Be  recoils  was  only 

^-35%.  Beryllium  metal  radiators  should  emit  detectable  fluences  of  recoils 

\2 

of  approximately  three  times  that  of  C  recoil  fluences  emitted  by  graphite. 
This  greater  sensitivity  of  beryllium  is  due  to  three  factors.  First,  the 

Q  ip 

elastic  scattering  cross  section  is  somewhat  greater  for  Be  than  for  C. 
Second,  the  lighter  Be  nucleus  carries  away  a  greater  fraction  of  the 
energy  from  a  collision  with  a  neutron.  Third,  the  smaller  the  Z  of  a 
particle  of  given  energy,  ■.".he  greater  is  its  range  in  matter.  Thus,  the 
sensitive  thickness  of  the  radiator  for  recoil  emission  is  greater  for 
lighter  particles.  However,  measurements  showed  that  these  advantages 
of  a  berylliwn  radiator  were  lost  due  to  the  insensitivity  of  the  Lexan 
track  recorder. 

12  16 

Experiments  showed  that  C  and  0  recoil  nuclei  were  efficiently 
recorded  by  Lexan.  Even  here,  however,  accuracy  of  measurement  was  less 
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than  desired.  The  etched  surface  of  Lexan,  although  acceptable  for  many 
applications,  did  not  prove  satisfactory  for  counting  of  short  recoil  tracks. 

A  substantial  and  somewhat  nonuni form  background  of  etch  pits  resulted  from 
small  surface  imperfections. 

A  survey  of  commercial  materials  did  not  reveal  a  recorder  with  the 
various  track  recording  characteristics  required.  It  was  subsequently 
experimentally  found  that  some  track  recorders  have  inherent  limitations 
in  their  ability  to  discriminate  between  ions  of  different  Z.  With  careful 
measurements  of  ions  near  the  registration  threshold,  it  was  found  that, 
in  Lexan,  Vy/Vg  increased  very  gradually  with  increasing  dE/dX.  In 
cellulose  nitrate,  Vy/Vg  increased  very  rapidly  with  dE/dX  once  the  registration 

threshold  was  reached.  It  is  this  second  type  of  track  recorder  which  is 

4  9 

required  to  discriminate  against  He  particles  while  still  recording  Be 

recoil  nuclei  with  maximum  efficiency. 

Cellulose  nitrate  plastics  of  the  usual  (commercial)  compositions  are 
4 

quite  sensitive  to  He  particles  of  energies  up  to  ^3  MeV.  However,  it  was 
found  that  the  sensitivity  may  be  varied  by  altering  the  kind  and  amount  of 
plasticizers  included.^  Plastics  low  in  plasticizer  content  tend  to  have 
poor,  nonuni form  etching  properties.  Also,  the  more  volatile  plasticizers 
tend  to  be  lost  from  the  surface  if  included  above  a  certain  level,  during 
preparation  of  plastic  sheets.  A  series  of  plastics  of  decreasing  sensitivities  was 
manufactured  by  loading  them  with  successively  greater  amounts  of  relatively 
nonvolatile  plasticizers.  By  irradiating  the  specimens  with  stopping  ^He 
particles,  it  wa*-  possible  to  select  plastics  whose  sensitivities  were  just 
below  that  needed  to  record  He  particles.  By  varying  fabrication  procedures, 
it  was  possible  to  arrive  at  a  plastic  recorder  which  also  exhibited  the 
desired  surface  characteristics. 


An  irradiation  with  14  MeV  neutrons  was  made  to  test  the  final  material. 
It  was  found  that  9Be  recoil  nuclei  were  registered  efficiently.  Very 
short  tracks  were  clearly  developed  and  longer  tracks  etched  out  with 
very  small  cone  angles.  Also,  the  surface  imperfections  did  net  appear  to 
be  a  hindrance  to  track  counting. 

A  series  of  beryllium  radiator-gold  degrader  units  were  then  manufac¬ 
tured  by  a  vacuum  evaporation  technique.  Aluminum  degraders  were  found  not 
suitable.  They  were  found  to  produce  background  tracks  even  in  plastics 
which  are  insensitive  to  ^He  particles.  These  are  very  short  recoi !  particle 
tracks.  The  track  densities  produced  were  considerably  smaller  than  from 
any  of  the  radiators  used.  However,  they  were  larq,:  enough  to  raise  the 
background  level  by  a  factor  of  ^2.5,  Where  thick  radiators  were  desired, 
gold  was  evaporated  onto  0.010  in  thick  beryllium  squares.  A  1  cm  square 
contained  four  different  gold  thicknesses.  For  thin  radiators,  a  tantalum 
sheet  of  0.010  in  thickness  was  used  as  a  backing.  Beryllium  films  of 
0.250  mg/cm^  were  evaporated  on  1  cm  squares.  Gold  was  then  evaporated 
onto  the  beryllium  as  before.  Table  1  gives  a  list  of  the  detector 
combinations.  For  neutron  energies  up  to  20  MeV,  any  beryllium  radiator 
over  1.9  mg/cm  is  thick  in  the  sense  that  recoils  of  the  maximum  possible 
energy  have  ranges  less  than  this  value. 


III.  Detector  Response 

A  method  has  been  developed  for  the  calculation  of  recoil  nuclei 
track  densities  per  unit  incident  neutron  fluence  for  the  HEND  system. 

The  physical  parameters  of  any  given  detector  can  be  inputted  to  a  computer 
program,  FLUENZ,  which  yields  response  versus  energy  for  a  discrete  grid 
of  neutron  energies.  It  is  assumed  that  neutrons  are  incident  normally 
on  the  face  of  the  detectors. 

The  geometry  of  the  system  is  shown  in  Figure  1  where  the  degrader 
is  omitted.  Since  the  neutrons  are  incident  normally,  the  recoil  fluence 
is  symmetrical  with  respect  to  the  scattering  angle,  e.  An  annular 
volume  increment  at  distance  s  from  P  is  given  by 

dV  =  2n  s2  sine  ds  de.  (1) 

The  number  of  nuclei  in  the  volume  increment  is  then  given  by 

dN  =  6.025  x  1023  xpxf  ,  (2) 

where  p  is  the  density  of  the  radiator  material.  If  the  radiator  is  not 
a  single  element  and  isotope,  each  nuclear  component  which  contributes  to 
recoil  fluence  must  be  calculated  separately.  The  proper  fractional 
density  of  the  isotope  is  substituted  into  Equation  2  for  p.  The  total 
elastic  cross  sectional  area,  in  cm^,  of  nuclei  within  the  volume  increment 
for  a  monoenergetic  neutron  beam  is  then  given  by 

dz  =  dN  x  a  x  1(T24 


dz  =  0.6025p  o  ~  y 


(3) 


where  a,  the  elastic  scattering  cross  section,  is  given  in  bams.  There¬ 
fore,  the  number  of  elastic  collisions  in  the  volume  increment  per  incident 
neutron/cm  is 


dV 


dC  =  0.6025p  o 


(4) 
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when  the  total  attenuation  of  the  neutron  beam  in  the  radiator  is  assumed 
to  be  negligible. 

When  an  elastic  collision  takes  place,  the  probability  of  a  recoil 
emission  per  steradian  at  angle  f  is  assumed  to  be  proportional  to  cose 
at  all  neutron  energies.  In  fact,  there  is  some  structure  in  the  angular 
scattering  probabilities  for  discrete  neutron  energies.  However,  these 
calculations  were  made  for  the  average  response  of  a  detector  for  0.5  MeV 
wide  energy  intervals.  For  the  neutron  energy  region  up  to  20  MeV,  the 
angular  structure  was  expected  to  average  out. 

The  recoil  fluence  from  a  single  collision  and  in  the  direction  of 
s  is  then  given  by 

F<;  =  £5®  (5) 

s  ns2 

The  fluence  in  the  direction  of  the  incident  neutrons  is  then 

Fn  =  Fs  cos  e 

F  .  co£e  (6) 

n  ns2 


Finally,  the  differential  response  at  P  from  the  volume  increment  is 

dR  =  dC  x  Fn 

Jn  0.6025p  a  dV  „  cos2  e 
dR  *  - S - x  Ts3- 


Substituting  for  dV 

dR  =  1.205p  a  sine  cos2  ede  ds/A 
The  analytical  solution  for  the  total  response  is  then 


R  =  V.?P$P--£  f7  f  sin  e  cos2  e  de  ds. 


o  Jo 


, 


After  the  first  integration  we  have 


^  _  1.205p  p 


s__w  sin  ecos2e  de. 


(7) 


[»! 


(8) 


The  quantity  s;flax  is  a  complex  function  of  e.  Physically,  it  is  the  maximum 
thickness  of  radiator  material  which  a  recoil  particle,  emitted  at  an 
angle,  e,  to  the  incident  neutron  beam,  can  penetrate  to  produce  a  detectable 
track  in  the  track  recorder. 


It  can  be  seen  that  s^  is  dependent  upon  initial  recoil  energy, 

and  thus  also  upon  scattering  angle,  e,  through  the  equation 

j.  _  4A  En  cos2  e 
r  ‘  (1+A)2 


Er» 

(9) 


where  En  is  the  neutron  beam  energy  and  A  is  the  atomic  number  of  the 

recoil  nuclei .  It  is  further  dependent  upon  the  following: 

R(E,Zn),  the  range-energy  relationship  of  the  recoil  ions  in  the 
radiator, 

tR,  radiator  thickness, 

R(E,7d),  the  range-energy  relationship  of  the  recoil  ions  in  the 
degrader, 

tD,  degrader  thickness,  and 

a E  ,  the  effective  minimum  energy  for  production  of  detectable 
tracks  in  the  track  recorder. 

For  the  neutron  energies  and  recoil  nuclei  in  question,  it  may  be  assumed 
that  no  upper  limit  on  recoil  energy  eyists  for  track  registration. 

However,  in  order  to  be  detected,  a  track  must  have  some  minimum  length. 
This  factor  is  responsible  for  the  minimum  energy,  AEr. 

The  substitution  of  an  analytical  expression  for  smax  in  Equation  8, 
for  complet  on  of  the  integration,  is  obviously  impractical.  Instead,  a 
numerical  integration  was  done.  In  order  for  the  numerical  integration  to 
closely  approximate  the  analytical  case,  the  integrand  had  to  be  calculated 
at  intervals  between  the  limits  of  integration  which  correspond  to  small 
changes  in  value.  The  value  of  s^x  is  at  a  maximum  at  e  =  0  and  decreases 
smoothly  to  zero  at  some  angle  less  than  =  n/2.  Therefore  calculations 
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at  1  degree  intervals  were  considered  sufficient. 
Equation  8  can  be  replaced  by 

r  .  rw  sine  cos2e. 

0=0 


(10) 


The  quantity,  K,  is  a  factor  necessary  to  renormalize  the  summation  to  the 
integral  value.  Since  the  original  limits  of  integration  were  0  to  n/2, 
then  for  90  intervals 

K  =  (n/2  -  0)/90 
K  =  0.01745 

Equation  10  then  becomes 

.........  qn 

01) 


R  =  0.021Q3p  a  £  -n  Q  2  6 

*  e=o  mx 


Since  range-energy  relationships  are  usually  given  in  terms  of  mass  thick¬ 
ness,  mg/cm^,  it  is  convenient  to  use  these  units.  The  factor  10" ^/p 
converts  mg/cn 2  to  an.  Thus 


2.103  x  10“5  a 

R  =  - J - 0Iosmax  sinecos2e» 


(12) 


and  s„v  is  now  to  be  calculated  in  units  of  mg/cm^.  Equation  12  holds 

moA 

only  if  the  radiator  is  composed  of  a  single  element  with  a  single  isotope. 
Otherwise  the  two  p‘s  are  not  identical  and 

R - rr2 - -  0^osmax  sin  6  cos2  6  (12*b> 

where  px  is  the  density  of  the  target  nuclei  within  the  radiator  and  p2 
is  the  total  density  of  the  radiator. 

The  sumnation  and  evaluation  of  Equation  12  was  done  with  Program 
FLUENZ  at  a  series  of  discrete  neutron  energies  and  for  selected  radiator 
materials  and  radiator  and  gold  degrader  thicknesses.  The  quantities 
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inputted  to  the  program  were  A,  atomic  number  of  the  scattered  nuclei,  a 
range-energy  table  for  the  radiator  material  and  gold,  called  R(I),  E(I), 
RA{ I ) ,  the  radiator  and  degrader  thicknesses,  called  TR(I)  ar.d  TD(I),  the 
grid  of  neutron  energies,  called  EF(I),  and  the  target  nuclei  elastic 
scattering  cross  sections  corresponding  to  the  neutron  energies,  called 
SIGMA(I).  Also  a  quantity  called  DELT  was  read  in.  This  was  a  thickness 
of  gold  equal  to  the  range  of  a  recoil  ion  of  energy  AEf,  the  minimum 
energy  for  production  of  a  track  in  the  plastic.  By  inputting  DELT  instead 
of  AEf,  an  extra  range-energy  interpolation  was  avoided,  thus  saving 
machine  time. 

In  Program  FLUENZ  a  nest  of  DO  loops  insured  that  the  variables,  EF, 
TR,  TD,  and  e  are  encountered  in  their  proper  order.  First  tables  of 
cose,  called  B(N),  and  cos2  6,  called  BTU(N),  were  calculated.  Then,  the 
elastically  scattered  recoil  energy,  called  EE(J),  is  calculated  according 
to  Equation  9.  Next,  the  recoil  ion  energy  necessary  to  penetrate  the 
degrader  and  register  in  the  plastic  is  found.  This  quantity  is  inter¬ 
polated  from  the  range-energy  table.  It  is  the  energy  corresponding  to 
gold  thickness,  WAX,  where 

DMAX  =  TD/B(J)  +  DELT 

and  B(J)  is  cose.  In  the  subroutine,  the  interpolated  quantity  is  always 
called  RANGE,  whether  it  is,  in  fact,  a  range  or  an  energy.  Next,  ranges 
in  the  radiator  corresponding  to  recoil  ions  of  two  different  energies 
are  interpolated.  One  energy  is  that  previously  found,  the  minimum 
energy  for  penetrating  the  degrader  and  registering  in  the  plastic.  The 
other  energy  is  EE(J).  The  second  range  minus  the  first  is  called  BIT. 
Therefore,  BIT  is  equal  to  the  maximum  range  in  the  radiator  which  a  recoil 
ion  can  span  and  still  penetrate  the  degrader  and  register  in  the  plastic. 
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The  maximum  thickness  of  the  radiator,  RMAX,  is  then  found  from 


RMAX 

= 

TR/B(0). 

If 

RMAX 

> 

BIT 

then 

SMAX 

* 

BIT. 

If 

RMAX 

< 

BIT 

then 

SKAX 

= 

RMAX. 

Thus,  Sj^x  is  equal  to  the  maximian  possible  recoil  ion  range  for  track 
registration  unless  this  range  is  greater  than  the  radiator  thickness. 

In  that  case  Sj^  is  equal  to  radiator  thickness.  In  the  program 
is  called  S.  It  is  transferred  into  the  subscripted  variable,  PATH(J), 
then  multiplied  by  sin  eccs2e  and  summed  over  the  90  angular  intervals  as 
in  Equation  12.  After  the  summation.  Equation  12  is  completed  by  a 
multiplication  by  the  proper  factor.  The  product  is  stored  for  printout 
as  FLU(I,K,L) .  The  subscripts  correspond  to  the  values  selected  for  EF(I), 
TR(K) ,  and  TD(L) . 

For  EF(I) ,  a  grid  of  40  neutron  energies,  separated  by  0.5  MeV  and 
extending  to  20  MeV,  was  employed.  The  corresponding  elastic  scattering 
cross  sections  used  were  averaged  over  0,5  MeV  intervals  about  the  neutron 
energies.  The  calculated  detector  responses  should  the«^efore  be  considered 
histograms  summed  over  the  same  energy  intervals  rather  than  true  differen¬ 
tial  responses.  The  cross  sections  of  9Be,  12C  and  160  collected  for  this 
purpose^*®  are  given  in  Table  2.  Some  gaps  occurred  in  available  data, 
especially  above  15  MeV,  and  extrapolations  or  interpolations  were  made  to 
fill  in  the  regions.  The  values  of  TR(K)  and  TD(L)  were  those  of  actual 
detectors  such  as  are  given  in  Table  1.  The  value  of  &Er  was  taken  to  be 
0.36  MeV.  This  value  is  difficult  to  determine  exactly  from  observations 
of  tracks  due  to  recoils  of  mixed  energies.  It  will  be  fixed,  finally. 
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TABLE  2 

Interval -Average  (0.5  MeV)  Elastic 
Scattering  Cross  Sections  Used  In 
Program  FLUENZ.  The  Units  Are  bams. 


Neutron 

Energy 

(MeV) 

9Be 

0.5 

4.2 

1.0 

3.6 

1.5 

2.26 

2.0 

1.85 

2.5 

2.23 

3.0 

2.87 

3.5 

1.76 

4.0 

1.38 

4.5 

1.22 

5.0 

1.165 

5.5 

6.0 

6.5 

^■11 

7.0 

1  1.027 

7.5 

8.0 

8.5 

!  0.995 

9.0 

9.5 

0.985 

10.0 

0.980 

10.5 

0.976 

11.0 

0.971 

11.5 

0.966 

12.0 

0.961 

12.5 

0.956 

13.0 

0.952 

13.5 

0.947 

14.0 

0.942 

14.5 

15.0 

15.5 

0.928 

16.0 

0.923 

16.5 

0.918 

17.0 

0.913 

17.5 

18.0 

«■■■!»?  til 

18.5 

0.899 

19.0 

0.894 

19.5 

0.889 

20.0 

0.835 

2.20 

3.15 

1.71 

1.82 

1.60 

0.880 

1.873 

1.488 

2.138 

2.995 

2.000 

2.350 

1.689 

1.537 

1.288 

1.249 

1.044 

1.426 

0.939 

1.400 

1.081 

1.449 

[■WMgSf 

1.043 

0.837 

0.689 

0.786  j 

0.575 

0.477 

0.695 

0.534 

0.755 

0.584 

0.628 

0.823 

0.665 

0.818 

0.700 

0.812 

0.733 

0.804 

0.748 

0.797 

0.744 

0.788 

0.738 

0.729 

0.773 

0.722 

0.766 

0.715 

0.759 

0.709 

0.752 

0.703 

0.745 

0.698 

0.739 

0.693 

0.733 

0.688 

0.727 

0.683 

0.722 

0.67 « 

0.717 

0.675 

0.712 
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by  conparing  calculations  of  detector  response  with  experimental  measure-* 
ment.  The  range-energy  tables  for  ions  in  radiator  and  degrader  materials 
were  taken  from  Northcliffe  and  Schilling.^ 

Figures  2  and  3  show  the  results  of  calculations  for  the  beryllium 
radiators  and  gold  degraders  given  in  Table  1.  Rather  than  present 
histograms,  smooth  curves  have  been  drawn  through  neutron  energy  grid 
points.  A  listing  of  the  program  is  given  in  Appendix  A. 
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IV.  Spectral  Unfolding 

The  development  of  a  computerized  unfolding  program  to  obtain  neutron  spectra 
from  HENO  measurements.  Program  HENQU,  is  in  progress.  The  unfolding  employs  an 
iterative  technique  similar  to  the  SAND  Program  of  McElroy,  et.al.®  The  SAND  unfolds 
activation  foil  measurements  into  neutron  spectra.  It  employs  a  sat  of  activation 
foils  with  overlapping  responses  (interaction  cross  sections)  covering  the  desired 
neutron  energy  range. 

An  inspection  of  the  problem  reveals  the  basis  for  the  unfolding.  Any  discrete 
neutron  spectrum  will  be  sampled  by  each  activation  foil,  according  to  the  cross' 
section,  to  yield  a  set  of  measured  activations.  A  model  spectrum  is  then  folded 
into  each  cross  section  curve  to  yield  a  corresponding  set  of  calculated  activations. 
A  comparison  of  the  calculated  and  the  measured  activations  gives  information  on 
how  the  model  spectrum  should  be  changed  to  approximate  the  experimental  spectrum. 

In  practice  an  iterative  scheme  generates  successive  approximations  to  the  experi¬ 
mental  neutron  spectrum  until  a  "best  fit"  is  found  for  the  activation  measurements. 
If  there  were  no  significant  errors  present  in  either  the  cross  sections  or  experi¬ 
mental  activation  measurements,  and  if  the  neutron  energy  range  were  well  covered 
by  overlapping  activation  foil  cross  section  peaks,  then  this  method  would  yield 
excellent  results.  In  practice,  fairly  good  results  are  obtained  over  wide  energy 
regions  (±10  to  ±30%  by  McElroy). 

The  least  reliable  of  McElroy's  measurements  tend  to  be  at  the  higher  neutron 
energies.  There  are  few  well  placed  interactions  of  high  cross  section  with  pro¬ 
per  decay  times  in  this  region.  Therefore,  the  HEND  system  should  3'ield  better 
results  within  its  region  or  applicability  (-^5  to  20  MeV)  than  do  activation  foils. 

A  summary  of  Program  HENDU  follows.  The  series  of  detector  response  curves 
may  be  thought  of  as  a  matrix,  R(I,J),  where  I  denotes  the  detector  grid  and  0  the 
neutron  energy  grid.  Each  detector  yields  a  measurement  of  track  density.  A, 
and  the  series  of  track  densities  may  be  denoted  A(I).  The  neutron  energy  grid 
consists  of  intervals  0.5  MeV  wide.  The  interval  center  energies  are  denoted 


Another  array  quantity  necessary  as  input  to  the  unfolding  program  is 
the  first  approximation  neutro:*.  spectrum.  The  first  approximation  should 
be  as  close  to  the  actual  experimental  spectrum  as  possible.  The  solutions 
to  iterative  fitting  programs  such  as  HEhlDU  are  non-unique  and,  in  general, 
the  less  "work"  necessary  to  gain  a  fit  to  the  experimental  spectrum,  the 
better  that  fit  will  be.  This  first  approximation,  CALC(-.'),  which  is  either- 
calculated  or  selected  by  other  Keans,  will  initially  be  inputted  to  the 
program  ’n  tabular  form. 

The  program  can  be  outlined  with  the  following  quantities: 

A(I) 

E(0) 

R(I,J) 

CALC(vJ) 

IMAX,  the  total  number  of  detectors  used, 

JMAX,  the  total  number  of  energy  intervals,  and 

KMAX,  the  maximum  number  of  iterations  allowed. 

The  iterative  portion  of  the  program  begins  with 

JMAX 

AC(I )=l  CAtC(J)*R(I.J)  (13) 

J=1 

where  AC(1)  are  the  calculated  track  densities  corresponding  to  the  measured 
track  densities  A(I).  In  fact,  if  CALC(O)  is  equal  to  the  actual  experimental 
neutron  spectrum,  then  within  the  statistics  of  the  measurements,  AC(I)-=A(/). 

CALC(J)  and  AC(I)  are  then  normalizee  to  the  measurements.  That  is 


IMAX 

CALC(J)  =  CALC(O)*  l  A{I)/[AC(I)*IMAX] 
1=1 


04) 


AC(I,J)  =  CALC(J)*R(I,J),  and 


05) 


JMAX 

AC( I )  =  l  AC(I,J)  (16) 

J=1 


Since  these  quantities  are  recalculated  over  and  over  again  by  the 
iterative  cycle,  they  can  be  ranked  according  to  iteration  number,  K,  for 
better  understanding  of  the  procedure.  Thus,  the  initial  approximation 
spectrum  may  be  labeled  CALC(J)^  and  calculated  track  densities  AC(I)^ 
since  no  iterations  have  yet  occurred. 

Next  the  quantities  DEV(I)^  and  STDV^  are  calculated.  DEV(I) 
are  the  individual  deviations  of  the  ratio  of  measured  to  calculated  track 
densities  for  each  detector  from  the  average  ratio,  in  absolute  fraction. 
This  is 

°EV(1)  -  (#&)  /  «T  #fj))  /  MMJ  -1-  07) 

STDV  is  the  current  iterative  value  of  the  standard  deviation  of  all  the 
DEV(I).  Thus 

riel  imax  rn 

STDVlkj  a  SQRT  (  l  DEV(I)L  J  **2  /  (IMAX-1)).  (18) 

1=1 

The  weighting  factors,  HATE  (J),  are  then  calculated.  The  first 
approximation  spectrum  is  transformed  into  the  first  iterative  spectrum  by 


CALC(J)^  =  CALC(0)£°J  *  EXP  (WATE(J)^). 


for  the  Re¬ 


iteration 


WATE(J) 


where: 


(19) 


(20) 
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(21) 


W(I,J)[k:i  =  AC(I,J)[k] 

AC(I)^ 

In  the  calculation  of  WATE(J),  the  values  of  A(I)  cannot  be  allowed  to 

equal  zero,  since  In  0  =  -  Ambiguous  values  of  DEV(I)  and  STDV  would 

also  be  obtained  for  this  condition.  Experimentally  these  values  can  equal 

-4 

zero.  For  these  calculations,  an  arbitrarily  small  value,  1  x  10  ,  is 
substituted  for  any  zeros. 

CALC(J)^  is  then  transferred  back  into  Equation  1  and  the  second 
iteration  begins.  The  iteration  number,  K,  is  set  up  to  keep  a  running 
score  of  how  many  iterations  have  taken  place.  When  KMAX  is  reached,  the 
program  is  terminated.  Another  method  of  terminating  the  program  is  by 
comparisons  of  STDV  for  successive  iterations.  A  knowledge  of  counting 
statistics  allows  one  to  calculate  the  statistical  errors  associated  with 
A(I).  From  these  an  expected  standard  deviation,  XSTD,  can  be  calculated 
and  inputted  to  the  program.  When  STDV  =  XSTD,  then  a  fit  has  been  found 
to  the  experimental  neutron  spectrum.  In  either  case,  a  listing  is  made 
of  DEV(I)  and  STDV  for  every  iteration.  This  allows  any  serious  errors  in 
measurements  to  be  recognized  and  a  record  ot  convergence  to  be  examined. 

Program  HENDU,  as  developed  to  date,  is  given  in  Appendix  B.  Some 
additional  quantities  have  been  added  for  practical  reasons.  A  table  of 
definitions  precedes  to  program. 

Up  to  the  present  time  Program  HENDU  has  been  tested  with  a  response 
matrix  R  (I,J),  composed  of  several  response  curves  of  the  thick  radiator 
type.  Simple  model  spectra  have  been  unfolded.  The  unfolding  took  place 
in  a  generally  satisfactory  manner.  With  few  detectors  having  widely 
spaced  neutron  energy  thresholds,  and  beginning  with  a  flat  first  approx¬ 
imation  spectrum,  the  unfolded  spectrum  teads  to  fit  the  model  spectrum 
in  segments  corresponding  to  the  inter threshold  gaps.  Figure  4  demonstrates 
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one  such  unfolding.  The  arrows  at  the  bottom  denote  the  first  values  of 
each  of  the  seven  response  curves.  The  first  response  begins  at  5  MeV. 
The  unfolded  spectrum  can  extend  down  in  energy  only  this  far.  The  high 
energy  part  of  the  spectrum,  which  was  well  sampled  by  the  responses, 
converged  rather  rapidly.  At  low  energies,  where  sampling  was  meager, 
convergence  was  much  slower.  One  method  of  improving  the  situation  is  to 
add  detectors  whose  responses  are  not  peaked  so  heavily  toward  higher 
energies.  Thin  radiator  responses,  such  as  those  in  Figure  3,  have  this 
property.  Also,  the  thresholds,  instead  of  being  rather  evenly  spaced, 
can  be  shifted  toward  lower  energies.  The  optimum  number  of  detectors 
will  be  decided  from  tests  under  experimental  ci rcunstances  and  including 
statistical  counting  errors. 
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V.  Applications 

The  HEbD  system  can  be  applied  to  the  spectral  measurement  of  high 

energy  neutrons  over  a  rather  wide  range  of  fluences.  The  saturation 

fluence  with  beryllium  radiators  is  about  10^  ng/cm2,  although  detectors 

with  the  more  sensitive  responses  may  begin  to  saturate  at  smaller  values. 

This  range  can  be  extended  by  using  thin  polyethylene  radiators.  The 

1  ? 

smaller  probability  of  emitting  C  recoils  from  these  radiators  increases 

n  7 

the  saturation  fluence  to  about  5  x  10  ng/cm.  At  this  fluence,  background 
track  densities,  due  to  recoil  events  occurring  within  the  track  recorder, 
become  very  significant.  It  is  necessary  to  count  large  numbers  of  tracks 
to  reuuce  statistical  errors  in  the  measurements.  Consequently,  manual 
counting  becomes  too  laborious  to  undertake. 

The  lowest  detectable  fluence  depends  on  both  the  energy  spectrum  of 
the  neutrons  and  the  particular  detector  configuration  in  question.  For 
the  more  sensitive  detectors  with  beryllium  radiators,  the  average  sensi¬ 
tivity  is  about  1.5  x  10-8  tracks/neutron,  as  seen  from  Figure  2.  A 
fluence  of  1.5  x  108  n0/cm2  would  yield  102  tracks/cm2.  This  is  a  low 
density  for  counting  recoil  tracks.  At  present,  even  the  best  of  plastic 
track  recorders  contain  background  pits  whose  density  is  a  substantial 
fraction  of  this  number.  In  addition,  manual  counting  of  such  densities 
is  time  consuming  if  statistical  errors  are  to  be  kept  low. 

The  effective  neutron  energy  range  of  the  detector  is  about  5  to 
20  MeV.  Below  5  MeV  efficiencies  become  quite  low  for  the  most  sensitive 
detector  configurations.  Neutron  energies  above  20  MeV  will  certainly 
produce  tracks,  but  the  response  in  chis  region  can  no  longer  be  described 


by  a  simole  elastic  scattering  procedure. 

Methods  for  automatic  counting  of  tracks  have  been  used^8 


or  are 


being  developed11*^2.  For  the  most  part  these  were  developed  for  counting 
of  fission  fragment  tracks,  but  the  same  methods  can  also  be  applied  to 

counting  of  recoil  tracks.  The  adaptation  of  an  automatic  counting 
technique  would  greatly  facilitate  the  use  of  these  detectors  at  all  fluence 
ranges . 

To  summarize,  some  of  the  advantages  of  the  HEND  system  are:  (1) 
small  size  (VI  cm^)  and  light  weight  (a  few  grams),  (2)  stability  (track 
registration  characteristics  of  nuclear  track  detectors  remain  constant 
over  long  periods  of  time),  (3)  permanently  recorded  tracks  (can  be  read 
any  time  following  exposure),  (4)  insensitivity  to  X-  or  gamma-ray  back¬ 
grounds  below  MO6  rads,  (5)  absense  of  flux  effects  (very  high  intensity 
pulses  of  neutrons  can  be  measured  without  detector  response  fluctuations), 
(6)  independence  tf  reaction  half-life  constraints,  and  (7)  wide  neutron 
energy  range  coverage  (cross  sections  exist  for  threshold  measurements  of 
neutrons  with  energies  between  ^5  and  20  MeV).  The  detectors  are  passive 
devices  but  active  readout  could  probably  be  attained. 


VI.  Conclusions 

At  this  tii?e,  the  configuration  of  the  detector  has  been  established 
and  the  type  of  information  to  be  gained  from  their  use  is  well  understood. 
The  ability  of  the  system  to  yield  useful  spectral  information  is  indicated. 

A  lengthy  but  necessary  development  of  a  srecial  track  recorder  for 
efficient,  selective  registration  of  recoil  particles  has  been  completed. 
This  was  achieved  by  reducing  the  sensitivity  of  a  cellulose  nitrate 
plastic  below  the  registration  threshold  of  ^He  particles.  Calibration 
measurements  can  now  proceed  with  the  final  version  of  the  detectors. 

A  calculational  method  of  determining  neutron  energy  response  of  the 
detectors  has  been  completed.  Experimental  verification  at  a  few  energies 


will  validate  the  calculations. 

A  spectral  unfolding  program  has  been  developed.  It  has  been  shown 
to  yield  model  spectra  for  situations  similar  to  the  experimental  conditions. 
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Programs  FLUENZ  and  HENDU  were  written  in  Fortran  IV.  They  were  run 
on  a  Univac  Series  70  (formerly  RCA  Spectra  70)  computer.  This  machine  is 
not  comparable  in  speed  to  those  generally  encountered  for  scientific 
application.  The  programs  ran  for  about  280  and  170  seconds,  respectively. 


A.  Program  FLUENZ 

Main  program:  reads  and  writes  all  input  data  deck  cards  and  writes 
output,  controls  program  flow;  calculates  recoil  particle  fluence  leading 
to  production  of  tracks  in  nuclear  track  detector  for  a  given  set  of 
detector  parameters. 

Subroutine  NTERP:  performs  interpolations  on  range-energy  tables  of 
energetic  ions  in  selected  media  (selects  either  range  from  energy  or  energy 
from  range). 

The  program  variables  are: 

A,  atomic  number  of  the  recoil, 

OELT,  a  range  in  gold  corresponding  to  the  minimum  energy 

of  the  recoil  nuclei  for  track  registration, 

R(N),  E(N),  RA(N),  range-energy  table  for  the  recoil  nuclei  in  the  radiator 

and  gold  degrader, 

EF(I),  the  neutron  energies  (MeV). 

SIGMA(I)  elastic  scattering  cross  sections 

TD(<J),  degrader  mass  thicknesses  (mg/cm2), 

TR(K),  radiator  mass  thicknesses  (mg/cm2), 

B(90),  cose  at  1  degree  intervals, 

BTU(90),  cos2e  at  1  degree  intervals, 

EE{90),  recoil  nuclei  energies  versus  o 

FLU(I,K,L),  calculated  track-producing  recoil  fluences, 

PATH (90),  intermediate  step  for  summing  FLU  over  e. 

Program  FLUENZ 

DIMENSION  au  (40,12,2),  EE(90),  PATH(90),  8(90),  BTU(90) 

DIMENSION  R(23),  E(23),  RA(23),  TR(2),  TD(12),  EF(40), 
SIGMA(40) 

READ  (5,11)  A,DELT 
F0RMAT(2n0.3) 
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1 

2 

3 

4 

5  11 


6 

READ  (5,12)  ((R(I),E(I),RA(I)),I=1,23) 

7 

12 

FORMAT  (3F10.4) 

8 

READ  (5,15)  (TR(I) ,1=1 ,2) 

9 

READ  (5,15)(TD(I),I=1 ,8) 

10 

READ  (5,15)(EF(I),  1=1,40) 

11 

READ(5,15)(SI6MA(I),I=1 ,40) 

12 

15 

FORMAT  (8F10.3) 

13 

WRITE  (6,18)  A,DELT 

14 

18 

FORMAT  (1H0.2F6.3) 

15 

WRITE  (6,17)  ((I,R(I),£(I),RA(I)),I=1,23 

15 

17 

FORMAT  (1H  ,15,3F8.4) 

17 

WRITE  (6,19)  (TR(I),I=1,2) 

18 

WRITE  (6,19)  TD(I),I=1 ,8) 

19 

19 

FORMAT  (1H  ,F7.3) 

20 

WRITE  (6,16)  ((I,EF(I).SIGMA(I)),I=1 ,40) 

21 

16 

FORMAT  (IH  ,1G,2F7.3) 

22 

00  5  N=l,90 

23 

B(N)=C0S((N-1.J*  0174533 

24 

5 

BTU(N)=B(N)*B(N) 

25 

DO  210  1=1 >40 

26 

X=4.*A*EF(I)/((1 ,+A)*(1.tA)) 

27 

DO  210  X=1 ,8 

28 

DO  210  L=2,2 

29 

DO  200  J=2,90 

30 

EE(J)=X*BPJ(J) 

31 

V 

EFFECT  OF  DEGRADER  THICKNESS 

32 

DMAX=(TD(K)/B(J))+DELr 

33 

IF(0MAX.GT.9.797)  GO  TO  121 

9.797  (MG/SQ  CH)  IS  MAX  VALUE  IN  INTERP  TABLE  FOR  BE  ON  AU 


34  C 


35 

50 

36 

C 

37 

C 

38- 

39 

40 

41 

42 

43 

44 

45 

C 

46 

401 

47 

48 

C 

49 

101 

50 

51 

102 

52 

103 

53 

54 

55 

56 

57 

200 

53 

C 

59 

121 

60 

61 

62 

202 

CALL  NTERP  (E,RA,DMAX,  RANGE) 

‘RANGE1  IS  NOW  THE  MIN.  ENERGY  ION  THAT  CAN  PENETRATE  DEGRADER 

AND  REGISTER  IN  PLASTIC 

TEST=RAHGE 

CALL  NTERP  (R,E,TE$T,RANGE) 

RANGE1=RANGE 

CALL  NTERP  (R,E,EE(J), RANGE) 

BIT=RANGE-RANGE1 
IF(BIT.GT.O.)  GO  TO  401 
GO  TO  121 

EFFECT  OF  RADIATOR  THICKNESS 
RMAX=TR(L)/B(J) 

IF  (RMAX-BIT)lOl ,101 ,  102 

S  IS  SET  TO  EITHER  RMAX  OR  BIT,  WHICHEVER  IS  SMALLER 
S  =  RMAX 
GO  10  103 


S=3IT 

CONTINUE 

PATH(J)=S 

PATHl-PATH(J) 

IF(PATH(J).LE.O. )  GO  TO  121 
PATH(J)=PATH(J)*SIGMA(I)*SQRT(1.-BTU(J))*BTU(J) 

CONTINUE 

SUMMATION  OVER  ALL  ANGLES 
JMAX=J-1 

IF(JMAX.LT.2)  GO  TO  203 
DO  202  M=2,JMAX 
PATH(M)=PATH(M)  +  PATH(M-l) 


{.  „  i»Vt l*l tzazl 


63 


GO  TO  209 


64  203 

65  209 

66  210 


68  21 


71  13 


74  22 


77  14 

78  211 


84  300 

85  301 


88  302 


90  303 


PATH(JMAX)=0. 

FLU(I,K,L)=2.33E-6*PATH(JMAX) 

CONTINUE 
WRITE  (6,21) 

FORMAT  (1H0) 

DO  211  1=2,2 
WRITE  (6,13)A,TR(L) 

FORMAT  {1H1  ,'MASS  HUMBER'  ,  F4.0,7*  ,  F6.3,'  MG/CM**2') 
DO  211  K=1 ,8 
WRITE  (6,22) 

FORMAT  (1H0, 1  DEGRADER  ENERGY  FLUENCE'/ 

DO  211  1=1,40 

WRITE  (6,14)  L,K,I,TD(K),EF(I),FlU(I,K,L) 

FORMAT  (1H  ,  315,2F8.3,E14.6) 

CONTINUE 

STOP 

END 

SUBROUTINE  NTERP  (DIST,ENER,EION, RANGE) 

DIMENSION  DIST(23),ENER(23) 

J=1 

IF  (EI0N-ENER(«.1))301 ,302,303 
RMAX=DIST(J) 

EMAX=ENER(J) 

GO  TO  304 
RANGE=DIST(0) 

GO  TO  309 
RMIN=DIST(J ) 
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91 

EMIN*ENER(J) 

92 

0=J+1 

93 

GO  TO  300 

94  304 

CONTINUE 

95 

IF( EMAX-EMIN) 305 ,305 ,306 

96  305 

RANGE=0. 

97 

GO  TO  309 

98  306 

RANGE=RMIN+(  ( (RMAX-RKIN)*(EION-EHIN)  J/(EMAX-EHIN) ) 

99  309 

CONTINUE 

100 

RETURN 

101 

END 
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8.  Program  HENDU 


Input: 


Card  No. 

Format 

Quantity 

Interpretation 

1 

515 

KTRL(L),  L=1 ,5 

KTRL(1)=0 

Given  values  of  4(1) 
are  read  in. 

KTSvL(l  )>0 

A(!)  are  computed  using 
a  known  spectrum. 

KTRL(2)*0 

tio  graph  is  plotted 

KTRL(2)=1 

Final  iteration  is 
plotted. 

KTRL\2),L>2 

No  assignment  has  been 
given  to  these  values 
yet  and  thus  they  can 
be  u:ed  as  control 
cards  in  the  future  use 
of  the  program. 

2 

3Ib,F10.3 

KMAX 

IMAX 

JNAX 

EMAX 

#  of  Iterations 

#  of  targets 

#  of  energy  points 

Maximum  energy,  i.e. ,EMAX=E(JMAX) 

in  general 

3 

9F8-3 

CALC(J},0=1 ,JMAX 

Value  o"  initial  approximation 
spectrum  at  energy  points  EO). 

& 

1615 

START( I), 1=1, IMAX 

First  J  value  for  which  R(I,«J)  is 
to  be  given  for  each  I.  It  is 
assumed  that  R(I,J)=0  for  all 

J<  START(I). 

5 

4(F10.5,10x 

R(I,d),J=1,JMAX 

For  each  I  value,  the  value  of  the 
response  curves  R(I,J)  for  each 
point  J  from  START ( I )  to  JMAX. 

6  a) 

m(i)7o 

9F8-5 

Value  of  theoretical  true 
spectrum  at  energy  values  E(J). 

b) 

KTRL(1 )=0 
9F8-3 

A(  I), 1=1, IMAX 

Track  densities  corresponding 
to  R(I,J). 
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Output: 

1.  Value  of  R(I,<J)  for  each  I  in  ascending  order  of  J.  No  title  given. 

2.  Values  of  A(I).  These  are  only  given  if  they  were  computed  from  theoretical 
curve,  i.e.,  if  KTRL(1 )>0  (title  given). 

3.  Iteration  number  starting  at  KMAX-5  (title  given), 
i)  Standard  deviation  (with  title). 

ii)  The  value  of  CALC(J)  for  each  J,  i.e.,  the  value  of  the  computed 
spectrum  at  each  J  value  (with  title). 

lii)  The  value  of  AC(I)  for  each  I,  i.e.,  the  iterated  values  of  the  track 
densities. 

4.  If  KTRL(2)  =  1 

For  each  value  of  J  a  point  is  plotted  giving  the  value  of  the  final 
iterated  spectrum.  A  smooth  curve  of  the  initial  theoretical  spectrum 
is  also  given. 

N.B.  as  yet  no  consideration  has  been  given  to  the  use  of  KTRL(1)=0  in 
which  case  the  smooth  curve  would  have  no  meaning.  Modifications  are 
1  needed  to  be  made  to  the  program  to  exclude  the  smooth  curve  when 
KTRLv » 


1 

2 
■> 

4 

5 

6 

7 

8 
9 

10 


PROGRAM  HENDU 
COMMON  A(20),R(20,40) 

COMMON  TH(42),CALC(42),E(42) 
DIMENSION  KTRL(5),AC(20) 

INTEGER  START (20) 

READ(5,16)  ( KTRL(L) ,L=1 ,5) 
READ(5,13)  KMAX,  IMAX,  JMAX,  EMAX 
READ(5,14)  CALC(J),  J=1,JMAX) 

READ  (5,72)  START  (K),  1=1, IMAX) 
DO  15  1=1, IMAX 


11  JMIN=START(I) 

12  JJ=JMIN-1 

13  DO  70  J=1,JJ 

14  70  R(I,J)  =  0.0 
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i 

i 

i 

i 


i 

\ 

i 


5  a 


I 

I 

& 

X 


4 


i  -3 


1 

$ 


a 


■M 


3 

3 


15 

15 

READ  (5,71)  (R(I,J),  J=JMIN,  JMAX) 

16 

DO  80  1=1 ,IMAX 

17 

DO  81  J=1 ,JMAX 

18 

81 

R(I,J)=R(I,J)*10**3 

19 

80 

CONTINUE 

20 

DO  73  1=1 ,IMAX 

21 

73 

WRITE(6,14)  (R(I,J),  J=1,JMAX) 

22 

IF  (KTRL(l))  18,18,  17 

23 

17 

READ  (5,57)  (TH(J),J=1 ,JMAX) 

24 

DO  31  1=1  ,IMAX 

25 

DO  32  J=1 , JMAX 

26 

IF  (R( I, J)-0. 0)32,32, 33 

27 

33 

A( I )=A( I )+TH( J )*R( I , J ) 

28 

32 

CONTINUE 

29 

31 

CONTINUE 

30 

GO  TO  19 

31 

18 

READ  (5,14)  (A(I),I=;,INAX) 

32 

19 

WRITE  (6,35) 

33 

WRITE  (6,12)  (A(I),I=1 ,IMAX) 

34 

DO  47  J=1 , JMAX 

35 

47 

E(J)=EMAX/UMX*J 

36 

C 

CALCULATION  OF  F 

37 

K=0 

38 

200 

F=0 

39 

DO  100  1=1,  IMAX 

40 

100 

AC(I)=0.0 

41 

DO  101  1=1 ,IMAX 

42 

DO  102  J=1 ,JMAX 
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43  102 

AC(I)=AC(I  )+CALC(J)*R(I,J) 

44 

IF(AC(!))  101,101,120 

45  120 

F-r+A(I)/AC(I) 

46  101 

CONTINUE 

47 

F=F/IMAX 

48  C 

NORMALIZING  INPUTS 

49 

DO  103  1=1 ,IMAX 

50  103 

AC(I)=AC(I)*F 

51 

DO  104  J=1 ,JMAX 

5?  104 

CALC(J)=CALC(J)*F 

53  C 

COMPUTE  STD  DEV 

54 

STD=0 

55 

DO  108  1=1,  IMAX 

56 

DEV=A(iy(AC(I)*F)-l 

57  108 

STD=STD+DEV**2 

58 

STDV=SQRT  (STD/{IMAX-1 ) ) 

59 

L=KMAX-5 

60 

IF  (L-K)  202,202,203 

61  202 

WRITE  (6,10)  K,STDV 

62 

WRITE(6,11) 

63 

WRITE(6,12)  (CALC(J),  J=1,JMAX) 

64 

WRITE(6,12)  (AC(I),  I='i,IMAX 

65  203 

K=K+1 

66 

IF(K-KMAX)  201,201,900 

67  C  CALCULATION  OF  WEIGHT  FACTORS 

68  201 

DO  106  J=l,  JKAX 

69 

DEN=0 

70 

TUM=0 
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72 

' 

IF(AC(I)-.00C1 )  107,107,109 

73 

109 

W=CALC(J)*R(I,J)/AC(I) 

74 

DEN=DEN*W 

75 

TUM=TUM+W*ALOG ( A ( I )/AC( I ) ) 

76 

107 

CONTINUE 

77 

IF  (DEN-.00001)  110,110,111 

78 

no 

WATE=0.G 

79 

GO  TO  106 

80 

111 

WATE=TUM/DEN 

81 

106 

CALC(«J)=CALC(J)*EXP(WATE) 

82 

GO  TO  200 

83 

10 

FQRMAT(/12H  ITER.  STEP= ,13,10X,5HSTDV=,F10.3) 

84 

11 

FORMAT (1 OX, 7HCALC(J)} 

85 

12 

FORMAT (  5F15.3) 

86 

13 

FORMAT (  I5,F10.3) 

87 

14 

FORMAT (  9F8.3) 

88  16 

F0RMAT(5I5) 

89 

71 

F0RMAT{4(F10.5,10X)) 

90 

72 

F0RMAT(10I5) 

91 

900 

IF(KTRL(2))901 ,902,901 

92 

901 

CALL  GRAPH 

93 

902 

STOP 

95 

END 

96 

SUBROUTINE  GRAPH 

97 

COMMON  IMAX,  JMAX,  A(20),R(’0,40) 
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9*1 

93 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 


COMMON  TH(42),CALC(42),E(42) 

DIMENSION  IBUF(IOOO) 

CALL  PL0TS(IBUF,1 000,8) 

CALL  PINIT('H200I ,'117', 20,‘ CALCULATED  SPECTRUM1) 

CALL  PL0T(0  0,0. 75, -3) 

I1=JMAX+1 

I2=JMAX+2 

CALL  SCALE(E,10.0,JMAX,1) 

CALL  SCALE(TH,5.0,JMAX,1) 

CALL  AXIS(0.0,0.0,1 3HENERGY  IN  MEV,-13,10.0,0.0,E(I1),E(I2)) 
CALL  AXIS(0. 0,0.0  £HSPECTRUM,8,5.0,90.0,TH(I1),TH(I2)) 
CALC(I1 )=TH(I1 ) 


CALC(I2)=TH(12) 

CALL.'LINEf  E,TH,UMAX,1 ,0,0) 
CALL  LINE(E,cALC,UMAX,L,-l,2) 
CALL  PL0T(1 5. 0,-0. 5,999) 
RETURN 
END 
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